A Cusp in QED at g = 2 
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We explore nonperturbative properties of the dimension-4 QED allowing a gyromagnetic ratio 
j ^ jd = 2. We determine the effective action V e a for an arbitrarily strong constant and homo- 
geneous field. Using the external field method, we find a cusp as a function of the gyromagnetic 
factor g in a) the QED bo-renormalization group coefficient, and in a b) subclass of light-light 
scattering coefficients obtained in the long wavelength limit expansion. We discuss precision QED 
results indicating an opportunity for resolution of known theory-experiment disagreements. We 
show possibility of asymptotic freedom in an Abelian theory for certain domains of g. 

PACS numbers: 12.20.-m, ll.15.Tk. 12.20.Ds, 13.40.-f 



Motivation: No known particle has exactly the Dirac 
value gD = 2 of the gyromagnetic ratio g. Determination 
of the higher order vacuum fluctuation correction to g 
provides the most precise test of Dirac-QED (D-QED), 
and this test is facing current challenges: 1) There is a 
three standard deviation discrepancy between theory and 
experiment involving the directly measured muon gyro- 
magnetic ratio P, Q , and theoretical result in 8th order 
a 4 [3[. 2) In a measurement of the muonic hydrogen 
Lamb shift, a discrepancy between D-QED and experi- 
ment is reinterpreted as a revised size of the proton [j, Q , 
inconsistent by five standard deviations with other avail- 
able experimental information [(| . 

In our opinion these recent experimental developments 
signal need to reexamine the framework of D-QED fo- 
cusing on the electron-magnetic field interaction. Recall 
that the D-QED pcrturbativc expansion is for the Dirac 
value g = gr> = 2. However, the theory of charged parti- 
cles interacting with the photon field, which we refer to 
generically as QED, taken as a stand-alone theory has a 
point-like electron or muon where in general g ^ 2 due to 
modifications by other interactions. Perturbative expan- 
sion around g — gjy is not appropriate should this value 
g = gu be a singular point. The aim of present work is to 
show that the there is a signularity at g = go, to study 
the nature of this singularity, and to lay foundation for 
a framework allowing exploration of \g\ > go- 

To achieve these goals we consider the extension to 
g 2 based on the renormalizable dimension-4 ac- 
tion 0, We find a singularity at g = go employ- 
ing the external field method: this means that we study 
the vacuum properties in presence of external constant 
and homogeneous electromagnetic fields, integrating out 
fluctuations of spin- 1/2 particles with g go- The re- 
sulting effective potential V e s is a generalization of the 
Heisenberg-Euler-Schwinger (HES) effective action [§-[l3j 
to arbitrary value of g. Our result is regular for all 

m<5d m 

For \g\ > go the HES effective action derived in proper 
time formulation [l4j based on dimension-4 renormaliz- 
able QED formulation becomes singular. We propose a 
natural extension for all values \g\ > go which shows 
that g = go (and other periodic recurrent values) is a 
cusp point as function of g. This extension resolves the 



known difficulties in the theoretical framework of g =/= go 
theories [l5|. Considering the beta- function and light- 
light scattering coefficients computed below, we discuss 
how the results we obtain for j / go can be tested by 
experiment and may impact at the required level the al- 
ready described two challenges QED faces today. 

Introducing magnetic moment |g| ^ go: One way to 

account for |<?| 7^ <7d is to complement the Dirac action 
with an incremental Pauli interaction term <5ju (a ■ B + 
id ■ E) = bp, <j a /3F a P/2 where E, B are the electromag- 
netic fields, F a P the electromagnetic field strength ten- 
sor, a a p = (i/2)[ r y a , 79] with 7 a the usual Dirac matrices, 
and a, and a = 7517 are the Pauli-Dirac matrices. This 
incremental Pauli interaction is a dimension 5 operator, 
[tjja a pF al3 ip] = L~ 5 . The coefficient 5fi consequently has 
dimension length, which in the case of a composite parti- 
cle such as the proton is naturally related to the particle 
size. This Dirac-Pauli (DP) equation has been a popular 
and effective tool to describe to lowest order the mag- 
netic moment dynamics of a composite particle of finite 
size, e.g. proton. 

A distinct approach to introduce an effective action for 
1.9 1 <?d is obtained by adding the full Pauli interaction 
term to the Klein- Gordon action 
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where n Q = id a + eA a . Note that the dimension of the if) 
field is [ip] = L^ 1 and consequently the Pauli interaction 
is dimension 4. We refer to the study of QED based on 
Eq. (TTJ) as g-QED, and the dynamical equation following 
from Eq. as the Klein-Gordon-Pauli (KGP) equation. 
g-QED is the s = 1/2 case in the study of particles of 
all spins in the Poincare group framework initiated by 
Rarita and Schwinger [l6j]. For recent developments see 
references in introduction to Ref. ■ 

Since there arc at least two distinct paths to introduce 
g / 2 corrections into relativistic particle dynamics, the 
question is in what sense these could be equivalent and 
if not, which of the two forms is appropriate for study 
of particle dynamics and/or vacuum structure and under 
what conditions: 

1) The DP approach, involving a dimension-5 operator, 
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requires new counter terms in each order. This is limiting 
DP approach to situations in which the physical particle 
properties are known and vacuum fluctuations need not 
be considered. Even so, vacuum fluctuations and the ef- 
fective action V c g have been considered in the DP i.e. 
modified D-QED approach (TJlll- 

2) In t/-QED the magnetic moment is point-like and con- 
sideration of g 7^ (?d does not require a higher dimen- 
sioned operator. Therefore the quantum field theory is 
renormalizable 0, H| , requires a finite number of counter 
terms, and vacuum fluctuations can be considered in any 
higher order. It should be remembered that g-QED con- 
stitutes an expansion around g = and not around 
g = ±2 as is the case for DP approach. Properties of 
the KGP-originating effective action were considered for 
general spin in Ref. [2(|, but this work did not recog- 
nize the restricted validity domain of the perturbativc 
approach, for spin- 1/2 — 2 < g < 2 . 

3) . A recent discussion of quantum field amplitudes with 
an anomalous moment [2l| also arrives at a second-order 
effective theory, but for a reduced two-component spinor. 
In view of the derivation and properties of their effective 
theory, a relation between KGP and DP can at best arise 
in an infinite order resummation in some specific appli- 
cations and not in general. 

Considering that Eq. ([T]) is 2nd order in time and has 
four components, the number of dynamical degrees of 
freedom present in Eq. (QJ is 8. That is, there are twice 
as many as in usual Dirac theory. For the case g = 2 
Eq. |T]) can be presented as the square of the operator 
j 5 D, D = 7 Q (i9 Q + eA Q )-TOand7 5 = i7°7 1 7 2 7 3 I 7| = 1 
is the 5th Dirac matrix. This means that for g = 2 
Eq. HI comprises exact duplication of the Dirac degrees 
of freedom. For g =/= 2 one must search for a projection 
restricting the full Hilbert space to the physical states. 

Vcltman has considered reduction of the number of 
dynamical components working in a two-component for- 
mulation. However, there are unresolved challenges 15 1 
in particular related to self-adjointness of the resulting 
spectrum and thus conservation of probability in tempo- 
ral evolution. By individually characterizing states, we 
will present another resolution of this problem that works 
in presence of externally applied fields and providing a 
result within the external field method. 

Eigenvalue-sum periodicity as a function of g: Our 

initial objective is to identify the physics content of the 
8 degrees of freedom and separate the Hilbert space into 
two equal size parts that each individually comprise a 
complete set of states at a fixed given value of g. To do so, 
we consider the Landau-orbit spectrum of the operator in 
brackets in Eq. JT]) in the presence of a constant magnetic 
field B 



E n = ±\Jm 2 +p 2 z + Q\eB\[(2n + l)Tg/2], Q = ±1, 

(2) 

where p z is the one dimensional continuous momentum 
eigenvalue and n is the Landau orbit quantum number. 



We have made explicit the presence of 8 eigenvalues for 
each value of B, corresponding to three different possi- 
ble choices of the signs. There are the usual two roots in 
Eq. (|2|) , a known feature of relativistic dynamics also seen 
in the Landau spectrum of the Dirac equation where the 
negative energy states become positive energy antiparti- 
cle 'hole' states. The =F factor inside the root in Eq. ^ 
arises from two possible particle spin projections onto 
magnetic field, corresponding to the spin degeneracy. 

Having recalled these usual features, we turn our at- 
tention to the new feature and write the eigen energy 
Eq. © in the form 



K = 



El 



\eB\ 



= Q[(2n+l)=Ffl/2], Q = ±l. (3) 



This exhibits a new spectrum duplication related to two 
possible values of Q. The quantity K is shown in the 
top portion of figure [T] as function of g. We see that 
between — 2 < g < 2 there is an exact duplication of the 
spectrum corresponding to Q = 1 and Q = —I. The 
'squared' Dirac operator produces two eigenstate-space 
copies which can be separated in particular applications. 
These are two sectors of the Hilbert space with the same 
physical content, and the Q = — 1 eigenvalues can be 
omitted. Thus for — 2 < g < 2 the effective action is 
obtained by the usual pro cedure, and the results have 
already been presented [I4( ■ 

For \g\ > 2, new physics content arises for exter- 
nal fields of any strength, including arbitrarily weak. 
First we note that taking Eq. ^ expression at face- value, 
naively some eigenstates could have E 2 < m 2 , which im- 
plies existence of bound localized states in the presence of 
a constant magnetic field. Such solutions are not required 
for completeness and would violate Lorcntz symmetry; 
for these reasons, such states cannot be admitted in the 



spectrum. This situation differs from the m 



Pz 







limit, in which states having K < signal Nielsen- Olesen 
instabilities of the conventional vacuum state [22[ ■ In the 
context of spin- 1 charged massless gluons in presence of 
color magnetic field, this instability has been used to ob- 
tain the leading coefficient of QCD rcnormalization group 
/3- function [23| . In g-QED for finite mass in the weak field 
limit, there is no magnetic instability. 

To compute the effective action we must define which 
states contribute to the physical spectral sum. The first 
step is to accomplish (like for the case \g\ < 2) separation 
of the Hilbert space into two sectors. We divide the states 
according to whether K > or K < and denote the 
respective sectors . The limit K = where two states 
coincide occurs at g = 2 since the KGP operator can 
be written as exact square of the Dirac operator. This 
situation recurs with the shift of g by 4fc, k £ Z. There is 
no change in the number of states in each of the Hilbert 
space sectors /C ± as an equal number of single particle 
states is exchanged between both sectors. 

The principle we use to determine which states en- 
ter the spectral sum is that there should be no localized 
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bound states in a constant magnetic field. In the no- 
tation just introduced, we require K > and the ]C + 
sector is chosen as representing the physical spectrum. 
This is an extension from the regular case \g\ < 2, where 
the usual procedure sums over the Q = +1 states and is 
equivalent to summing over the JC + state space. Seeing 
as K > implies E 2 > m 2 , the physics is a continuous 
extension of the case g = 2, for which it is proved that 
E 2 > m 2 for arbitrary magnetic fields, i.e. there are no 
bound states [24J. 

Looking far outside the principal domain — 2 < g < 2, 
we see that relativistic Landau eigenstates cross between 
/C ± at each = 2 + 4k, fceZ. As the graphic represen- 
tation top frame of Fig.[T] shows, for each of the Hilbcrt 
space sectors K.^ we have periodicity of the Landau lev- 
els a function of g. Therefore, the sum E n over K. + 

leading to the real part of V e s{B 2 ), is a periodic function 
of g, a result we will find explicitly. This periodicity does 
not apply to individual Landau eigenvalues as is seen in 
Eq. ([2]) . In computation of vacuum fluctuations the trun- 
cation of the Landau cigenstate n-sum to any finite value 
breaks the periodicity as well. 

The choice of JC + as the physical state space has clear 
advantages and resolves the challenges encountered by 
Veltman fl5j | : In addition to maintaining self-adjointness 
of the KGP system, it makes the quantum field theo- 
ries based on semi-spaces JC^ each individually unitary, 
because the number of states is conserved in transiting 
through the singular points e.g at \g\ = 2, and for \g\ > 2 
we omit the localized solutions. Moreover, our proposal 
makes the spectrum and by extension the quantum the- 
ory a continuous and analytic extension from the domain 
\g\ < 2. Our approach preserves translation invariance 
of the vacuum, which would be broken by any localized 
bound states in the constant-ficld-fillcd vacuum. It is 
critical to note that had we separated the sectors along 
the sign of Q, the contents of the theory would be dif- 
ferent for \g\ > 2 and unitarity would be violated since 
the 'wrong' levels would be included in the physical half- 
space. 

Effective action for |g| < 2: We briefly summarize re- 
sults for \g\ < 2 [l4j|, as these are needed to understand 
the present case of \g\ > 2. For constant fields the effec- 
tive action is manifestly covariant, and can be written 
as a function of the Lorentz-invariant field-like quantities 
a, b 

b 2 -a 2 = B 2 -E 2 = ^F af) F afs = 25, (4) 

(ab) 2 = (E ■ Bf = Qi^c^A^) 2 = V 2 , (5) 

where ±a are electric-field-likc and ±ib arc the magnctic- 
ficld-likc eigenvalues of F Q/9 . a is considered electric- like 
because a — > \E\ on taking the limit b — ¥ 0, and similarly 
b -> \B\ in the limit a — s- 0. 

The Schwinger-Fock proper time method to evalu- 
ate the effective action exploits properties of the 'squared' 
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FIG. 1: Top: Squared eigenvalues Eq. @ of KGP in magnetic 
field; the solid and (blue) dashed lines are for Q = +1, and 
respectively (— ) and (+) spin eigenvalue; the dotted and dash- 
dotted (red) lines are for Q = — 1, and respectively (+) and 
(— ) spin eigenvalue. Bottom: coefficient functions: /i.o(ff) as 
defined in Eq. (|I2[) and fo,2 and f^,o as defined in Eq. (|I3[) . 
Two full periods are shown. The values of g where the sign 
of the functions fij changes is indicated. 

Dirac equation and thus it can be used to study arbitrary 
value of g. The effective action can be written in form 

V eS - ^ ^ e-^ m2 -^F(eau, ebu, §). (6) 

of Jo u " 

For g = 0, 2 the proper time integrand F(eau, ebu, g) 
was reviewed in Ref . [IJ| . The generalization throughout 
the interval \g\ < 2 is accomplished by inserting into 
Schwinger's Eq. (2.33) in the last term a co-factor g/2 
leading to 

2 sinha; siny 2 

The subtraction —I in Eq. removes the field- 
independent constant. The logarithmically divergent 
charge renormalization term is isolated and discussed be- 
low. Note that Eq. (|6]) would be divergent for \g\ > 2 if 
Eq. (|7|) were to be used in this domain. 

Effective action for |g| > 2: To extend Eq. to 
\g\ > 2, we consider in more detail the eigenvalue summa- 
tion method we introduced above, follo wing the work of 
Heisenberg and Euler H and Weisskopf [Tflf. The math- 
ematical tool used was the L. Euler summation formula, 
leading to the Bernoulli functions B2k(x) and Bernoulli 
numbers &2k = ^2fe(0). The sum of the Landau energies 
Eq. © involves the form J2 n f( x + n )- L. Euler developed 
the technique for such sums, which manifest an integer 
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shift symmetry in the variable x — > x+n' [25, 26]. Due to 
this shift symmetry, the Bernoulli functions B 2 k(x) that 
arise in the context of L. Euler summation of Landau en- 
ergies E„ , Eq. (|2|) are the periodic Bernoulli functions, 
given by the Fourier series |27| 

R (+\-( Uk-l ( 2fc ) ! y cos(27rnt) 
B2k(t)-{-l) z^-T 1^ (n7r \2k > ( 8 ) 

n=l V / 

(here only needed for an even value of index, 2k). In 
the unit interval, < t < 1, the periodic Bernoulli 
functions are equal to the Bernoulli polynomials, e.g. 
B 2 {t) = B 2 (t) = t 2 - t + 1/6, < t < 1. Outside 



The new mathematical element on the RHS in Eq. © 
is that to assure the necessary periodicity we introduced 
in accordance with Eq. ((5J) a series of Bernoulli functions 
with t = g/4 + 1/2. We note that the right hand side of 
Eq. (JSJ) agrees exactly with the known expansion [29| in 
the domain of Eq. (|7j! \g\ < 2. This expression provides 
an analytical continuation into the domain \g\ > 2 having 
the periodicity property of the effective action identified 
in study of the full set of eigenvalues. 

Each term in Eq. © produces a well-defined result for 
all g upon performing the proper time integral Eq. ([5]) . 
The form Eq. Q is thus a unique and convergent exten- 
sion to \g\ > 2 determined by the Euler summation of 
the eigenvalues Eq. <j2j) . Even after the removal of the 
charge renormalization subtraction term (first term on 
RHS of Eq. @) the remainder of the effective action is 
manifestly periodic in g but not in e. 

Nonperturbative in g renormalization group f3 
function: The first non constant term on the right hand 
side of Eq. © proportional to a 2 — b 2 isolates the loga- 
rithmically divergent one- loop 0(a) V e g subtraction re- 
quired for charge renormalization. The coefficient of this 
term is related to the /3-function coefficient bo as is dis- 
cussed e.g. in section 5.1 in Ref. The two next 
terms (a 2 — b 2 ) 2 and (ab) 2 correspond to lowest order 
effective field-field interaction potentials describing light- 
light scattering. Setting in the remainder denominator 
on RHS of Eq. (|10[) y = produces next term in the 



the unit interval, the periodic Bernoulli functions Eq. ([5]) 
B2k{t) repeat the polynomials' behavior on < t < 1 in 
each subsequent period. 

Dividing the Landau energies by 2\eB\ to make the 
coefficient of n unity, we see that t — > g/4 + 1/2 and 
hence we recognize that the periodic Bernoulli functions 
with argument t = g/4 + 1/2 appears in the effective 
action, arising from the summation of eigenvalues. The 
explicit representation of the argument of Eq. ([5]) in terms 
of Bernoulli functions is arrived at employing the analytic 
transformation of the integrand of Eq. © [28|, [2j| . 



(9) 



(10) 



expansion, etc. 

We now consider explicitly the running of the coupling 
constant a within the g-QED loop expansion of the /3- 
renormalization function 

»->%• «-)-^ + ^ + .... ("J 

The first sum in Eq. ©, for g = 2, £^=1 l/(7m) 2 = 1/6 
and implies the value of bo = —4/3, where factor 4 in- 
dicates the 4 components of spin-1/2 particle. For arbi- 
trary g, bo(g) is obtained using Eq. ([5]) to identify this 
sum as B<x{g/4 + 1/2). The character of this function is 
manifest by reconnecting periodic domains of the familiar 
Bernoulli polynomial B 2 (t ) = t 2 — t + 1/6 and the result- 
ing bo(g) coefficient is given in each domain g £ [fffc-i, gk] 

6o = -|A,oOr) = -| (§(» - 4*) a - 1) , (12) 

where /i,o(<?) is shown in bottom frame of Fig.[TJ The 
subscripts of fij indicate the powers of the Lorentz in- 
variants in polynomial expansion ji,jS l V^ in Eq. ((9|). We 
sec in Fig.[T]that as a function of g, the Dirac value gu = 
±2 is an upper cusp point with /i.o(<?) < /i,o(2) = 1- 
For clarity, two periods are shown in Fig. [T] 

Note that our result arises from the KGP equation 
applying a nonperturbative method in g to one loop ex- 
pansion. This approach is necessary in order to obtain 
the behavior of the /3-function for \g\ > 2. At g = ±2 we 
find the unexpected cusp. This feature is missing in per- 
turbativc consideration of /3(g) at one loop level which 



F( X ,y°) = (x 2 -y 2 )2± 1} + 2 £) + 1} 



n=l 



(nn) 2 ^— ; (nir) 4 



\n=l V ' \n=l 



\ 2 \ 

cosn-7r(§ + l)\ \ q 



(tut) 2 J I ' - v^' 1 " 2 
where we separated the lowest powers of fields from an exact remainder function 
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produces the same functional dependence on g as seen 
in Eq. (fT2)) setting k — 0. As our study shows, a pertur- 
bative expansion around g = has a finite convergence 
interval \g\ < 2. 

The following implications for g-QED of the properties 
of the renormalization group coefficient &o(s) are note- 
worthy - we address in the following discussion the range 
of values of g shown in Fig.[TJ 

1. ) We recognize g is an independent 'large' coupling con- 
stant. The first-order D-QED expands around g = ±2, 
which points arc identified as being non-analytic in the 
<7-QED framework. 

2. ) For any value of g not at the cusp the magnitude 
|&o| decreases (and thus the speed of 'running' decreases) 
compared to its value at g = 2. Considering that the 
coefficient of the magnetic spin term in Eq. ([l} is dimen- 
sionless there is no new scale appearing in association 
with g. 

3. ) The presence of the cusp in bo implies that the run- 
ning coupling of g-QED, a(q) comprises the cusp as well. 

4. ) A cross check and confirmation of our result for 60(5) 
is obtained in perturbative domain considering the limit 
g — > where bo{g — > 0) differs only by a minus sign and 
the number of degrees of freedom from the known behav- 
ior of scalar particle 'QED'. The minus sign is due to the 
commutation relation needed in closing the fcrmion trace 
in loops, whereas it is absent in scalar boson loops. 

5. ) In the principal domain \g\ < 2 the functional depen- 
dence on g we find agrees with the result Eqs. (53-57) 
seen in Ref . [1| . Specifically, the leading term for large q 2 
of the vacuum polarization function, evaluated within the 
framework of g-QED is — abn(g)/(2ir) ln(— q 2 /m 2 ), seen 
explicitly in Eq. (55) of Ref. ||. 

6. ) As the above limit shows, for a range of appropri- 
ate gyromagnetic moment values g (including 5 = 0) 
bo(g) > is possible. This produces asymptotic freedom 
behavior for fcrmions interacting alone with an Abelian 
charge. The switch between the infrared stable and the 
asymptotically free behavior occurs in the principal g- 
domain twice, at g = ±2/vo = ±1.155 and continues 
periodically e.g. for g — 4 — 2/^/3 = 2.845. This mecha- 
nism of asymptotic freedom generation by g-driven sign 
reversal is implicit in Eq. (56) of Ref. Q (valid in prin- 
cipal domain \g\ < 2), but the new mechanism allowing 
Abelian confinement has not been recognized there. The 
values of g where the sign of the functions foj changes is 
indicated in Fig.[JJ up to periodic recurrence. 

Light-light scattering as function of g: We find that 
the cusp at \g\ = 2 reappears in a directly observable 
phenomenon inherent in the Heisenberg-Euler action, the 
light by light scattering. For the general case of both 
electric and magnetic fields present, using Eq. ^ we find 
up to fourth order in the fields 
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where both ^0 and /o,2 arc normalized to g = 2 val- 
ues and presented in Fig.[TJ Jo, 2 includes a product of 
two Bernoulli functions with cusp and so has a steeper 
cusp. Importantly /o,2 enters the V 2 term which one ac- 
tually measures in laser light scattering off a magnetic 
field [13, Hl|. In general, our finding is that all fij(g) 
for j > have cusps at g = 2 whereas all fi,o(g),i > 1 
are continuous and diffcrentiable at g = 2, being propor- 
tional to higher order > 2 Bernoulli functions that have 
vanishing derivatives at g = 2. Thus only coefficients of 
terms involving powers of the pseudo scalar field invari- 
ant V 2 — (E ■ B) 2 display cusps at g = 2. 

Discussion: We found new physics arising for \g\ > 2 
for arbitrarily weak fields in (/-QED. We proposed a new 
cigenstate sorting based on sign of K, Eq. (0, leading to 
a self-adjoint theory that retains Poincare symmetry and 
contains a complete set of particle- antiparticle states, and 
thus preserves probability in time evolution and analyt- 
icity as function of g, up to a countable set of singular 
points. 

While Eq. ([7]) is an analytic function of g, the integral 
of Eq. Q with the proper time weight Eq. © does not 
exist for \g\ > 2. Thus a naive extension of HES effective 
action to \g\ > 2 is not possible. This parallels the ob- 
servation that the Klcin-Gordon-Pauli operator Eq. ([JJ is 
not self-adjoint for \g\ > 2. We have presented a careful 
study of how the eigenstate level crossing can be recog- 
nized and states assigned to half-spaces of the full Hilbert 
space, leading to a natural self-adjoint extension and a 
valid theoretical g-QED framework for \g\ > 2. The cusp 
and related nonperturbative in g effects arise from im- 
plementation of the self-adjoint extension described. The 
origin of the cusp is in the periodic crossing of eigenener- 
gies in the spectrum of Landau eigenstates seen in upper 
section of Fig. [1] showing the quantity K, Eq. (|3|). 

The top frame of Fig. Q] illustrates the case of a (weak) 
magnetic field only, a similar result is obtained for the 
case of an electric field only, leading to the unique form 
Eq. © . Within this expression we have shown cusps at 
fj = .9d for two physical quantities computed for arbi- 
trary g: 

• The renormalization group coefficient bg proportional 
to function /x,o, see Fig.[TJ 

• The light-by-light scattering in the long-wavelength 
limit comprising a smooth function /2,o ; an d for the term 
(E ■ B) 2 the cusp function /0.2, sec Fig. [TJ 
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These nonpcrturbative in g one-loop results require in- 
finitely many contributions from eigenvalues Eq. ([2]) , and 
for this reason the effects predicted here are only visible in 
phenomena arising from vacuum fluctuations, such as the 
/3-function and effective light-light interactions. We have 
checked that these results can be arrived at directly by 
the method of (^-function regularization following Weis- 
skopf 0. 

Our results agree in the fundamental domain —2 < 
g < 2 with earlier perturbative work: the functional de- 
pendence on g is explicit and the same for the vacuum 
polarization as had been obtained in Ref. [§| in Eq. (56). 
We have shown by explicit computation that an expan- 
sion around g = is valid for \g\ < 2 only. 

Conclusions and outlook: Difficulties of D-QED as a 
stand-alone theory have been known for some time, be- 
ginning with the work of G. Kallen [32j], and perturbative- 
D-QED is believed by many to be semi-convcrgent only. 
Exploration of g ^ 2 in a renormalizable theory requires 
the dimension-4 g-QED based on KGP equation. How- 
ever, g-QED has to begin with 8 degrees of freedom and 
appropriate division into two half-Hilbcrt spaces is re- 
quired. Restriction to the usual Dirac-likc 4 degrees of 
freedom is difficult, as a theory with g ^ 2 is in general 
not unitary [Tfl . We resolved this problem, and from the 
solution we discovered that the Dirac value g = gu = 2 is 
a cusp point of the effective action V e g, Eq. ([6]) evaluated 
in renormalizable g-QED approach. 

This finding implies that the D-QED expansion around 
g = gu could be incomplete at sufficiently high order. To 
see the problem, imagine that we partially resum g — 2 
diagrams with Dyson-Schwinger method finding an ef- 
fective electron with g > 2. In the next step we want to 
compute the vacuum polarization inserts in other g — 2 
diagrams. Attempts in D-QED framework will encounter 
new divergences as the g — 2 correction is dimension-5 op- 



erator. On the other hand, we can accomplish this task in 
g-QED: we use the non-perturbative in g renormalization 
group coefficient bo to characterize the vacuum polariza- 
tion loop insert and there are no new divergences. How- 
ever, the result contains the cusp, and thus is different 
from the finite order perturbative expansion of D-QED. 

We believe that the higher order vacuum polarization 
modification we described would be most visible for the 
Lamb shift of muonic hy drogen which is dominated by 
the vacuum polarization [33l435| . Similarly, we expect 
that our light-light scattering cusp modifies the corre- 
sponding contribution to the muon g — 2. We discussed 
in seperate work how the top-quark loop modifies the two 
photon [3(| and two gluon [33] decay of the Higgs. At- 
tempts to evaluate these results using D-QED methods, 
that is a g — 2 non-renormalizable extension, would fail 
in evaluation of the loop integrals which in general are 
divergent in this case. 

Our study shows how a complete theory of a point-like 
fcrmion with \g\ > 2 can be constructed within g-QED in 
order to allow dynamical description of real world spin- 
1/2 particles. We have obtained the HES effective po- 
tential for an elementary particle with gyromagnetic ratio 
g 7^ 2 nonperturbativcly in g, see Eq. ^ and Eq. ©. We 
demonstrated a cusp as function of g at the Dirac value 
g = go = 2. We have shown how this cusp enters the 
(3- function and {E ■ B) 2n terms of light-light scattering. 
An interesting theoretical consequence is the possibility 
of asymptotic freedom in an Abelian theory with anoma- 
lous magnetic moment originating in the reversal in sign 
of the renormalization group coefficient bo for g in specific 
domains much different from g = 2. 

Acknowledgments: We thank I. Bialynicki-Birula, and 
R. Stora for careful reading and valuable comments. This 
work was supported by a grant from the US Department 
of Energy, DE-FG02-04ER41318. 



[1] G. W. Bennett et al. [Muon G-2 Collaboration], Phys. 
Rev. D 73, 072003 (2006) 

[2] G. Venanzoni, J. Phys. Conf. Ser. 349, 012008 (2012) 

[3] T. Kinoshita and M. Nio, Phys. Rev. D 73, 053007 (2006) 

[4] A. Antognini, et al., Science 339, 417 (2013). 

[5] R. Pohl, et al, Nature 466, 213 (2010). 

[6] P.J. Mohr, B.N. Taylor, and D. B. Newell, Rev. Mod. 
Phys. 84, 1527 (2012). 

[7] C. A. Vaquera-Araujo, M. Napsuciale, R. Angeles- 
Martinez and , JHEP 1301, 011 (2013) 

[8] R. Angeles-Martinez and M. Napsuciale, Phys. Rev. D 
85, 076004 (2012) 

[9] W. Heisenberg and H. Euler, Z. Phys. 98, 714 (1936). 
[10] V. Weisskopf, K.Dan. Vidensk.Selsk.Mat.Fys.Medd. 24, 
3, (1936). 

[11] J. S. Schwinger, Phys. Rev. 82, 664 (1951). 

[12] M. Reuter, M. G. Schmidt and C. Schubert, Annals Phys. 

259, 313 (1997) 
[13] G. V. Dunne, In *Shifman, M. (ed .) et al.: From f ields 

to strings, vol. 1* 445-522 and at: [hep-th/0406216 . 



[14] L. Labun and J. Rafelski, Phys. Rev. D 86, 041701(R) 
(2012). 

[15] M. J. G. Veltman, Acta Phys. Polon. B 29, 783 (1998) 
[16] W. Rarita and J. Schwinger, Phys. Rev. 60, 61 (1941); 

M. Napsuciale, M. Kirchbach and S. Rodriguez, Eur. 

Phys. J. A 29, 289 (2006) [hep-ph/0606"308] . 
[17] R. F. O'Connell Phys. Rev. 1761, 1433 (1968) 
[18] W. Dittrich, J. Phys. All, 1191 (1978) 
[19] P. M. Lavrov, J. Phys. A18, 3455 (1985). 
[20] S. I. Kruglov , Annals Phys. 293, 228 (2001) 

[hep-th/0110061] . 
[21] A. J. Larkoski and M. E. Peskin, Phys. Rev. D 83, 034012 

(2011) 

[22] H. B. Nielsen and P. Olesen, Nucl. Phys. B 160, 380 
(1979). 

[23] N. K. Nielsen, Am. J. Phys. 49 (1981) 1171; N.K. Nielsen 
and P. Olesen,, Nucl. Phys. B144, 376 (1978). 

[24] P. Gornicki, L. Mankiewicz and R. Zembowicz, J. Phys. 
A 20, 6593 (1987). 

[25] L. Euler, Memoires de l'Academie Imperiale des Sciences 



7 



de St.-Petersbourg 5, 45 (1815); [Translation available 

atarXiv:0705.0768 [math.HO]]. 
[26] T. M. Apostol, Amer. Math. Mon. 106, 409 (1999). 
[27] Qiu-Ming Luo, Math. Comput. 78, 2193 (2009). 
[28] B. Muller, W. Greiner and J. Rafelski, Phys. Lett. A 63, 

181 (1977). 

[29] Y.M. Cho and D.G. Pak,Phys. Rev. Lett. 86, 1947 (2001). 
[30] G. L. J. A. Rikken and C. Rizzo, Phys. Rev. A 63, 

012107 (2001) and references therein; and F. Bielsa,et 

al. I arXiv:091 1.45671 [physics.optics] . 
[31] S.-J. Chen, H.-H. Mei and W.-T. Ni, Mod. Phys. Lett. 

A 22, 2815 (2007 ); and H.-H. M ei, W.-T. Ni, S.-J. Chen 

and S.-s. Pan, arXiv:0911.4776 [physics. ins-det]. 
[32] G. Kallen, Consistency problems in quantum 



electrodynamics, CERN- "Yellow report"-57-43, 
DOL10.5170/CERN-1957-043; and Quantum Elec- 
trodynamics, Springer- Verlag, New York 1972, 233p, 
ISBN 0-387-05574-6 

[33] A. Antognini, F. Kottmann, F. Biraben, P. Indelicate, 
F. Nez, and R. Pohl, Ann. Phys. 331, 127145 (2013). 

[34] U. D. Jentschura, Annals Phys. 326, 500 (2011) 

[35] A. Czarnecki, U. D. Jentschura and K. Pachucki, Phys. 
Rev. Lett. 95, 180404 (2005) [Erratum-ibid. 95, 199903 
(2005)]. 

[36] L. Labun, and J. Rafelski, arXiv:12 d9T046l [hep-ph] . 

[37] L. Labun, and J. Rafelski. larXiv:1210.3150l [hep-ph]. 



